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ABSTRACT 


■Non-thermal emission which occurs in the cores of the 9.4 and 10*4nm COg 
bands on Mars has been recently identified as a natural atmospheric laser. The 
emission is believed to be excited by absorption of near-IR solar flux, followed 
by collisional transfer to the 00°1 state of COg. A numerical model based on 
this mechanism is developed which includes the solar pumping contributed by 
v^ZxlO^ near-IR COg lines. The NLTE radiative transfer problem is solved to 
obtain the 00°1 vibrational state population. Tnis model successfiilly 
reproduces the existing center-to-1 imb observations, although higher spatial 
resolution observations are needed for a definitive test. The model also 
predicts total fluxes which are close to the observed values, . Jhe strength of 
the emission is predicted to be closely related to the instantaneous near-IR 
solar heating rate. 

Calculation of the COg level populations in this model supports the 
identification of this emission as a natural atmospheric laser. Both Mars and 
Venus are predicted to exhibit level inversions; at large zenith distances 
several percent of the emergent intensity will be due to stimulated emission. 
While the stimulated ejnission present in these lines is not large enough to be 
of importance to mesospheric radiative equilibrium, it has other interesting 
consequences. The construction of large volume radiation-pumped lasers, which 
utilize COg planetary mesospheres as a gain medium, is theoretically possible. 


I. INTRODUCTION 


strong non-thennal emission in the 9.4 and 10.4wm bands of CO 2 was first 

reported for the atmospheres of Mars and Venus by Betz (1976) and Betz et al . 

(1976), using infrared heterodyne spectroscopy. The emission occurs in the 

cores of the C0« absorption lines and is only detectable at spectral resolving 
6 ^ 

powers of >10 • The emission is believed to be formed by absorption of 
near-infrared solar flux, followed by collisional transfer of the absorbed 
quanta to the 00°1 level of CO^. and radiative decay at 10pm. This process may 
be enhanced in the presence of water vapor (Johnson et al. 1976). Recently, 
Mumma et al . (1981) have demonstrated that in the case of Mars this emission 
arises from a population inversion, making it a natural laser.. ^ Theoretical 
pumping calculations have been given by Johnson et al . (1976), but these 
calculations accounted only approximately for the absorption by near-IR CO 2 
bands and neglected the effect of radiative transfer in the 4.3 urn band. The 
understanding of this emission and its implications is still not complete, and 
several interesting questions remain to be answered. Although the effect of 
emission in. these bands on the heating and cooling budgets of the Martian 
mesosphere is believed to be small (Ramanathan and Cess 1974), can it be a 
useful diagnostic tool in the study of mesospheric thermal balance in CO 2 
atmospheres? To what extent is the level inversion important in this context? 
The emission is believed to be optically thin, however, the assumption of 
optical thinness may not be valid for the long atmospheric path lengths 
encountered by photons emitted at large zenith angles. In this case stimulated 
emission may have a significant effect on the emergent intensity, and could have 
an appreciable impact on the mesospheric cooling rate. This possibility is best 
considered through theoretical modelling, since current observations do not have 
sufficient angular resolution to impose meaningful constraints on excess 
limb-brightening due to stimulated emission. 

This paper presents a more complete calculation of the pumping of the laser 
emission by including the solar pumping due to individual lines in all 
significant near-IR CO 2 bands. We solve the NLTE radiative transfer problem in 
the 4.3 urn band to obtain the population of the upper vibrational state. The 
possibility of significant laser amplification at large zenith distances is 
evaluated. 
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II. THEORETICAL CALCULATIONS 

a. The atmospheric models 

Prior to deriving and solving rate equations for the vibrational state 
populations, we must *adopt model atmospheres for the mesospheres of Mars and 
Venus. Two essentially ad hoc model atmospheres were adopted for the Martian 
mesosphere (Fig. 1). The Mars models were assumed to be isothermal in the 
mesosphere, with temperatures bracketing the range of plausibility based on 
Viking measurements (Seiff and Kirk 1977). The lower temperature (120K) agrees 
best with our own measurements of the thermal width of the emission lines 
(Deming et al . 1983). The Mars models were taken to extend fr.o^n the surface, at 

o 

6.3 mbars pressure, to 10 mbars, in 100 layers. The Venusian model mesosphere 
(Fig. 2) was taken from the Pioneer Venus "morningside model" (von Zahn et al . 
1980) and day probe measurements (Seiff et al . 1980). The Venus model was 
assumed to extend from 10" mbars to a lower boundary at 3 mbars, in 100 layers. 
The models thus have an altitude resolution of •/’0.2 scale heights. The 
calculated laser band emission is insensitve to the details of the model 
atmospheres at pressures less than *^10'° mbars, but the models were extended to 
higher altitudes in order to provide an adequate boundary for the 4. Sum 
radiative exchange calculation (described below). The laser band emission is 
also completely insensitive to the temperature structure of the models at 
pressures below »/*0.3 mbars. Since we will be concerned with the behavior of the 
laser emission at extreme zenith angles, we compute optical depths using a 
spherical atmospheric geometry. 

b . The Rate Equation for 00° 1 

In modelling the pumping of the laser emission we adopt the basic 
mechanism outlined by Betz (1976) and Johnson et al . (1976). These authors have 
discussed and eliminated alternative excitation processes. The mechanism for 
the excitation of the upper state of the laser emission, the 00°1 state, is 
absorption of incident solar radiation in the near-IR bands with subsequent V-V 
collisional transfer of the quanta to 00°1. Direct excitation of 00°1 by solar 
absorption at 4.3um was neglected by Johnson et al . (1976), but we include this 
direct pump in, our calculations. The pumping process is illustrated 
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schematically on Fig. 3. In order tc calculate the source function for the 
laser emission it is necessary to establish rate equations for both the 00*^1 
and (02^0)31,11 levels. In establishing these rate equations we make 

several assumptions. First, we assume that the relative populations of the 
rotational levels are at their LIE values throughout the models. Since only a 
few kinetic collisions are required to thermalize the rotational levels, this 
assumption is valid at the pressurer> considered here. Second, we assume that 
all quanta, excited by absorption of solar quanta in near-IR combination 
bands, are instantaneously relaxed by V-V transfer to the fundamental state. 
Third, we assume that collisional exchange of vibrational quanta from minority 
isotopic forms of COg, to the normal isotopic form, occurs instantaneously, and 
that the reverse processes never occurs. The temperature profile of the 
atmosphere is regarded as fixed, and we do not couple our vibrational state 
calculations to a radiative equilibrium temperature model. 

Let us denote the vibrational population of 00°1 by No (cm"*^); we sinsilarly 
denote the vibrational population of C(10°0), (02°0)]j by We choose to cast 
the calculations for the lower levels in terms of the I level since the majority 
of observations have been made in this band. For convenience, we denote 
C(10°0), (Q2°0)]j as 10°0’. 


We 


can write the rate equation for N 2 as: 


Ax j 'i’y -3;;, Me -^^1 

) = M, [Zuj ^ (,) 

L L * ^ 

. i i 

In this formulation '' represents the mean intensity at frequency v in a 

band connecting vibrational states i and j ; is the emission profile of the 

band, and includes the individual rotational line profiles. A,, and B.. 

1 i 3 

represent the Einstein A and B coefficients for the vibrational transitions; in 
this general notation is understood to be zero wnen state j is at higher 
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energy than state 1. .iil^j^indicates the velocity-averaged collisional rate (in 

sec~^) from vibrational state j to state k, when induced by collision 
with a molecule in state i, i.e.: 

In the most general case we would also formulate similar rate equations for 
all other vibrational states, as well as radiative transfer equations for the 
bands connecting these states. Simultaneous solution of these coupled equations 
with the appropriate boundary conditions (e.g. solar flux incident on the top of 
the atmosphere) would yield all of the relevant level populations. However in 
the cases of interest here (primarily the mesospheres of Mars and Venus) a 
considerable simplification is possible. Keeping only the most significant 
terms from (1) we have: 

- U) -A V, ^ da) (3 

+• (<] 

J } Ep 

These terms represent spontaneous emission in the laser bands (a), V-T 
collisional de-excitation of 00°1 to C(10°0) ,(02°0)]j (b), spontaneous 

emission to the ground state in the 4.3 ^m 00^1-Q0°0 band (c) , V-T collisional 
excitation of Q0°1 from the [(10°0) ,(02^0)]j jj levels (d), absorption from the 
ground state to 00°1 (e), and V-V collisrsnal transfer to 00^1 from higher-lying 
levels having one or more asymmetric stretching quanta excited (f). The factor 
v(j) accounts for the fact that de-excitation of overtones may produce more 
than one molecule in 00°1. We note that V-T collisional excitation and 
de-excitation of 00^1 has been assumed tc occur through the [(10°Q),(02°0)]I,II 
levels. Oe-activation of 00°1 in laboratory studies occurs to levels with 
bending and symmetric stretching vibrations excited. These levels, in turn. 


PA0P- iS 

QUAUTY 


6 


ORIGfWAL PASS fS 
OF POOf< QUALITY 

transfer vibrational quanta among themselves with large cross sections (Moore et 
. al. 1967, Finzi and Moore 1975). The level C(10°0) ,(02°0)]I ,II is formally 
talcen as the site of coHisional transfer from 00®1, although in practice this 
transfer probably occurs through higher-lying levels such as 

We note that (3) contains no term representing either absorption or 
stimulated emission in the lOum laser bands themselves* since these bands are 
expected to be optically thin. We will examine and justify this assumption 
later. 


We can rewrite the mean intensity in (3) as: 






^0 




nz^ (0) d CO + o 


/ 

AO 


(^3 


a) 


where the first term represents the solid angle average of directly incident 
solar photons, and represents photons created by emission from 00 1. We 
can also rewrite the write the most significant terms in the 

rate equation for N. as: 

w 

/^4 J <1 j y (a^ -f- (L) 




o 


(J) 
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This formulation accounts for direct absorption to level j (a), collisional 
excitation of level j from the 00^1 state (b), spontaneous emission from level j 
”to the ground state (c) and V-V collisional de-excitation of to od®l (d). 

w 

Since the levels are specified to contain one or more asymetric stretching 
quanta, will be large (Finzi and Moore 1975). We thus rewrite this rate 
equation as 



This formulation is equivalent to the assumption that all overtone quanta 
excited by solar absorption are instantaneously relaxed to 00°1. The mean 
intensities can be written as before: 

OJ 

In the case of the near-IR bands the first term, which represents the solid 
angle average of incident solar photons, is dominant over the self-emission of 
the atmosphere 
for 00°1 as: 


represented by • We can therefore write the rate equation 


- hi, hi, hi, 11),^ J 




where the last term now represents both direct solar excitation of 00°1 and V-V 
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transfer to 00°1 from overtone levels which are excited by near-IR absorption 
of incident sunlight. Under conditions of interest here we may take the ground 
state population 11^ as equal to tt»e total CO 2 number density. Equation (8) will 
be coupled to a similar equation for through the term In practice, 

however, this coupling is not important and (8) can be solved for N 2 /NQ by 
adopting an LIE value for N^. This possible because at high altitudes the 
solar absorption term on the right side of (8) is dominant, and the N^NqQj^ 2 ° 
term is negligible. At lower altitudes, where the becomes 

important, collisional processes dominate in the corresponding equation for N^. 
Similarly, we can include the under tiie assumption that the lOum 

bands are optically thin. At the lowest altitudes in the models, optical deptn 
effects at lOum become significant, but at these altitudes the term N 2 N^ a° 2 i 
dominates the left side of Eq. (8). 


The principal complication in the solution of (8) occurs through the term 
containing • When this term is evaluated at layer i in the atmospheric 

model it represents the mean intensity in the 4. Sum band wnich is produced by 
emission from 00°1 in a11 other layers of the model. At layer i we write 




R 




^^■TT J 


CO 



Where e-. 1s defined as the probability that a pnoton emitted in the band at 
level j in the atmosphere will be emitted in the solid angle du, and transmitted 
to level i without being absorbed in intennediate layers. aX- is the geometric 

J 

thickness of layer j. This integral equation formulation of the NLTE radiative 

transfer problem follows Dickinson (1972, 1976). Similar calculations have been 

performed by Kumer and James (1974). We calculate values of ^ i , j = 

1,100. This calculation is numerically tedious because each of these 
4 

100x100=10 terras involves numerical integration over u, integration over the 
^100 rotational lines in the band and integration over «/’10'^ frequencies in each 
line. Moreover, Voigt profiles must be used for the rotational lines since the 
emission arises in the Voigt transition region of the mesosphere. We affect 
integration over u by assuming azimuthal symmetry and integrating over u = cos e 
by 6-point gaussian quadrature, e being defined as the angle between the zenith 
and the direction of photon propagation. 
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Having calculated the (8) reduces to a system of 100 simultaneous 
linear equations, which can be solved by matrix inversion or by iteration. In 
^e 4 >resent instenc^ ^ have -solved -thi^ “system by matrix inversion. In 
addition to this exact solution, we also obtain an estimate for N 2 by using a 
**cool-to-space" approximation (Rodgers and Walshaw 1966, Dickinson 1972). This 
approximation removes the coupling between layers and allows us to solve each 
equation (8) for N 2 in terms of only local quantities. We implement this 
approximation in equation (8) by ignoring the term hlc^ox. 
replacing N 2 A 2 Q where <e^^> is the angle-averaged 

probability that a photon emitted into the upward hemisphere escapes to the top 
of the atmosphere without further absorption. This "escape-to-space" 
probability is calculated in the same way as the exchange terms. It 

(jj ® . 

is illustrated for the Venus model on Fig. 2. The number densities obtained for 
00°1 at the pressures of relevance mbars) do not differ significantly 

when the result for the cool-to-spaci approximation is compared to the exact 
radiative exchange calculation. 

... In solving (8) for N2 we adopt A2^ = 0.8 sec"^ (Murray et a1 . 1974) and A2Q 
» 441 sec”^ from the recent high precision measurements by Rogers et al . (198 ). 
We take = 3.64xlO“^^T (cm^^ sec“^) from Rosser et al . (1969) and Moore et 
al. (1967). 

b. The Rate Equation for 10° 0 

As long as the emission is optically thin, the population of its lower 
state has little effect on the emergent intensity. The fact that no emission is 
seen in the cores of the isotopic lines argues strongly for optical thinness at 
disK center, but this argument cannot be used at the extreme limb because the 
observations have insufficient spatial resolution. Moreover, Mumma et al . claim 
the levels to be inverted and we wish to calculate the lower level population in 
sufficient detail to support or deny this claim. 


OF POOR QUALITY 


The population of 10°0 is closely tied to the population of 01^0 through 
^olllsional V-V -exchange. Rate egua'ti'on? for these two levels can be written 
under the cool-to-space approximation as: 






3l 


L<e'V 


c 

'C ~ ~I3 


HXo3 = 


1^3 li -^ 3 . + 

* T<s'“> A,. 


OD 


Here denotes the population (cm‘^) of 01^0. The first term in (10) is a 
solar pumping term. It occurs because solar excitation of \>2 levels is assumed 
to be followed by fast (i.e. instantaneous) V-V reactions such as: 


If V-V exchange between 10°0 and 01.^0 were also instantaneous, i.e. fl °^2 
were very large, then solar excitation of overtone levels could be 
assumed to populate 01^0 directly (e.g. Dickinson 1972). However, we used log 
q°j^2*"13.07+1.368x 10”^T sec'^) from Bulthuis (1973). This rate, although 

large, is not fast enough to be considered as instantaneous, especially at the 
highest altitudes in our models. Other measurements, however, have given a 
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value of a°j ^2 magnitude larger than the Bulthuis rate (e.g. Rhodes 

et al . 1968). We have therefore repeated the model calculations using an order 
of magnitude greater value for This decreases the 10°0 population 

scfnewhat, leading to slightly more optical amplification in the laser bands, but 
has no other significant effect. In calculating the summation in Eq. (10), we 


assume that only solar excitation of [(10 1) ,(02 1)] 
significantly. 


I, II 


contr ibutes 


Since we argued that in practice can be derived independently of N^, we 
can use the values of N 2 obtained from (8) in the solution of (10) and (11). We 
adopt log Q°-q * -14.91+2.15x10"^! (cm^^ sec"^) from Winter (1963). We used 
Aj^ 2*2.9 sec’^ from the band strength estimates of Kostkowski and Kaplan (1957) 
and A-3n*l«44 sec"^ from Rogers et al . (1982). The angle-avera'ged 

3Q 13 

escape-to-space functions <c > and <e > were calculated by numerical 

13 

integration, as described above. An initial calculation of <e > was made by 

assuming an LIE population for 01^0. We then solved (10) and (11) 

13 

simultaneously for N, and N^. Values of <e > were then recomputed using the 

, A j 1 

new 01 ’0 populations. This process was repeated until the change in the 01 0 
population ms less than 0.1%. Since the 10°0 population does not affect 01^0 ' 
stroiigly^ liris iterative procedure converges rapidly. 

In calculating the solar absorption terms in (8) and (10), the band 
profiles, , were represented as sums of non-overlapping rotational line 

profiles. The rotational lines so included were taken from the 1973 version of 
the AFGRL line tape (McClatchey et al., 1973). These lines were grouped by line 
strength and a single rotational line profile was computed as^ representative of 
each group. This profile was then replicated into the profile as many 

times as number of lines it represented. All of the CO 2 lines on the tape were 
included, even isotopic lines. Many of these weaker lines give an insignificant 
contribution to the solar pumping. As noted earlier, excitation of quanta by 
absorption in isotopic lines was assumed to transfer instantaneously to 
The resultant profiles represent approximately 2x10 individual rotational lines 
in the region from 2000 to 9500 cm"^. Using this procedure the solar flux 
absorbed by a group of lines could be quickly computed as proportional to the 
flux absorbed by the representative profile. This is a good approximation at 
the pressures of relevarice here, because the typical line spacing is large 
compared to the line widths and the lines can be treated as non-overlapping. In 
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Implementing this procedure it is necessary to recognize that this collective 
representation of line absorption will be temperature-dependent because of 
differing lower state 'energies for lines in a particular strength class. The 
distributions of line strength were accordingly recomputed for each atmospheric 
model t although a single distribution was used at all depths in each model. In 
the case of the Venus model, which is not isothermal in the mesosphere, this 
distribution was computed at I85K. A Voigt function was used in computing the 
rotational profiles, and a pressure broadening coefficient of 0.07 cm"^ atm”^ 
was adopted for all lines. Values of the incident solar intensity appropriate 
to Mars and Venus were taken from Allen (1973). The incident solar flux for 
Mars was based on its mean distance, the eccentricity of its orbit causes an 
appreciable variation (±19%) in this number. Calculations wer^ made at a number 
of values of the solar zenith angle. The pumping line strength distributions 
were compiled over *>‘1500 cm'^ intervals, the exact value depending on the band 
spacings, and a single value for the incident solar intensity was used for each 
interval. Measured in photons cm“^ sec’^ sr’^^ne solar intensity does not vary 
rapidly with frequency over the near-IR region. The most important pumping 
bands in each interval are listed in Table I. 


Having calculated the populations of 00°1 from (8) and of 10°6 from (10) 
and (11), we also calculate the emergent intensity and optical depths in the 
core of a particular rotational line of the 00°l-[(10®0) ,(02°0)] j band, usually 
the 10.33Mm R(8) line. 


III. RESULTS 

a. The pumping rates 

Evaluation- of the solar absorption terms in (8) was done separately for the 
four frequency intervals of interest, and the resultant pumping rates into 00°1 
are shown for the 120K Mars model on Fig. 4. We confirm the claim of Johnson et 
al. (1976) that the 2.7um bands, (see Table I) in the 3000-4500 cm“^ region, are 
the dominant source of pumping. However, because of the large amount of 
saturation at the altitude of the emission peak (*/'70 km), even the very weak 
lines in the 6000-9500 cm"^ region contribute within a factor of 3 of the 2.7um 
region. We also find that, contrary to a statement by Johnson et al., direct 
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solar excitation of 00”l is significant. In fact, at the highest altitudes in 
the model this 2000-3000 cm’^ region pumping is the dominant means of exciting 


b . The 00°1 and 10°Q level populations 


The populations we calculate for the 00°1 and 10^0 levels are given on 
Figs. 5 and 6 for the Mars 120K and 170K models. The corresponding calculation 
for the Venus model is given on Fig. 7. In all of these models the 00°1 level 
becomes inverted with respect to the IQ°Q level. We therefore confirm the 
claim of Mumma et al . (1931), that level inversion exists In the emitting 
region. This inversion is ubiquitous at high altitudes, it appears even in 
models which do not include the full set of pumping lines. In the 120K Mars 
model the inversion extends down to altitudes below the emission peak; this 
occurs because of the much lower thermal excitation rate for 1P°0 at lower 
temperature. Thus we find that the lower altitude boundary of the inverted 
region depends on the kinetic temperature and on other characteristics of the 
model. At high altitudes a’ level inversion is unavoidable in all models since 
collisional excitation of 10°0 becomes unimportant, while direct solar 
excitation of 00°1 is significant. Other processes, such as solar pumping of 
10*^0* are not sufficient to perturb the inversion, and the ratio N(00°l)/N(10°0) 
. N^/Nj approachas j ^ 7 • 

c. The height dependence of the emission 

-3 -1 

Figs. 5, 6 and 7 show the total number of photons emitted (cm sec “) in 
the 9.4 and 10.4|im bands as In the 170K Mars model the emission peaks at 

an altitude of 75 km and drops to half intensity at 97 and 57 km. Thii behavior 
is essentially identical to the height dependence calculated by Betz (1976) and 
Johnson et al . (1976), who also used a 170K isothermal model atmosphere. 

However at an altitude of >/'130 km Johnson et al . show a secondary peak in the 
emission. Our preliminary models also showed this secondary peak, but it 
disappears when larger numbers of pumping lines are included. The Mars 120K 
model shows emission which peaks at 67 km, and falls to half intensity at 85 and 
51 km. The altitude of the emissicn peak is thus rather insensitive to the 
temperature of the model. In the case of Venus we find that the emission peaks 
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at 109 km and falls to half intensity at 120 and 101 km. Johnson et al . found 
the Venusian emission to peak at a^ilS km, but they used a wholly different model 
atmosphere. This comparison again shows that the altitude for peak emission is 
insensitive to the- details of the input model atmosphere. This is a useful 
result because it implies that the emission can be reliably associated with a 
fixed layer in the mesosphere. Repetitive observations of the frequency width 
of the emission cores would seem to be an excellent method for remote monitoring 
of temperature variability in the mesospheres of Mars and Venus. 

d. Limb brightening and optical depth . 

Having calculated the 00°1 and 10°0 level populations, we are now in a 
position to calculate optical depths in the lines and to evaluate the degree of 
amplification which occurs because of the inversion. We have done this 
calculation for the core of the 10.33um R(8) line because it is one of the 
strongest lines in the band, and because Oeming et al . (1983) have observed it 
extensively. Figs. 5, 6 and 7 show the total optical depth for this line when 
it is observed at the sub-solar point, with the line of sight perpendicuVar to 
the solar direction. This optical depth is shown as a function of the minimum 
altitude along the optical path. In this computation we neglect the small 
variation in solar zenith distance along the line of sight and we also neglect 
refraction. This calculation corresponds to the case when the extreme limb of 
the planet is observed, and the subsolar point occurs exactly at the limb. At 
the extreme limb the optical depth for R(8) is negative, indicating gain in the 
line due to stimulated emission. In the case of the Mars 120K model the maximum 
gain occurs when the line of sight reaches a minimum altitude of 66 km; in this 
instance the optical depth is -0.07. This means that each photon which 
traverses the entire line of sight has a 1 in 14 cnance of producing a: 
additional photon due to stimulated emission. Since an emergent photon on 
average traverses only half the total distance, 4% of the emergent photons are 
due to stimulated emission. When the minimum altitude of the line of sight 
decreases to 56 km the optical depth increases to zero(i), i.e. absorption and 
stimulated emission exactly cancel. At this "transparent point", the region of 
the line of sight which lies below 60 km produces absorption, since the 
inversion disappears below 60 km. However this absorption is exactly cancelled 
by stimulated emission which occurs for those regions on the line of sight which 
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The relative intensity emergent from the top of the atmosphere in the core 
of 10.33nm R(8) is shown as a function of angle in Fig. 8 for the 120K Mars 
model. In this figure the maximum intensity emitted at large zenith angles is 
more than 20 times the intensity emitted in tne zenith direction. This increase 
in intensity is almost entirely due to the increase in optical path length for 
this optically thin emission; the small additional increase at the limb caused 
by stimulated emission is shown by the inset to tne figure. We let be the 
direction between the angle of emission and the zenith, as measured from the top 
layer of the atmospheric model. We expect that emission from a very extended 
plane parallel slab will peak at cos e^*0. The emission shown , on Fig. 8 peaks 
at cosegsrO.2 because of the curvature of the atmosphere. The Venusian emission 
is predicted to have a similar center-to-1 imb dependence, but it peaks at a 
smaller value of cos (0.13) because the Venus atmosphere has a larger radius 
of curvature. 

The total relative flux (integrated over solid angle) emergent in the 9.4 
and .10.4^m bands, F (sg), is shown as a function of solar zenith distance (e^) 
for the Mars 120K model on Fig. 9. Because of the sphericity of the atmosphere 
the emergent flux is not zero when the solar zenith distance is 90 degrees. 

Since the emission' is optically thin, we expect that the relative flux in 
individual rotational lines would show the same dependence on solar zenith 
distance. 

Having calculated the effect of stimulated emission we can now evaluate its 
importance to the physics of the atmosphere, and we can justify its otnission 
from the rate equation for 00°1. It is clear from Fig. 8 that stimulaved 
emission will have a negligible effect on the solid angle average of tne 
emergent intensity and we can therefore neglect its impact on the radiative 
cooling rates in the mesosphere. It is also clear that its omission from (8) is 
justified because when averaged over solid angle it is insignificant in 
comparison to spontaneous emission. Inclusion of a stimulated emission term in 
(8) would only decrease the 00°1 population below the values we nave calculated, 
leading in turn to less stimulated emission. We thus conclude that laser 
amplification in the 9.4Mm and 10.4nm bands, while present in the mesospheres of 
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Mars and Venus, is of no Importance in the radiative equilibrium of the 
atmosphere. In another context, however, the laser amplification which occurs 
in ttie^e lines is startlingly Targe, and has interesting, if somewhat 
speculative, consequences (see sec. g). 


e. The total flux emitted 


The total flux emitted by the Mars models depends on the temperature of the 

model. The 120K model predicts a total flux of 14.4 ergs cm“ sec" and the 

-2 -1 

170K model predicts 20.6 ergs cm sec . The Venus model predicts a total 
emitted flux of 75.3 ergs cm"^ cm"^ sec”^. All of these values refer to the 
flux emitted with the Sun at the zenith, i.e. F^(0). The modeU predict that 
the total emergent flux should increase as the kinetic temperature increases. 

At higher temperatures the radiative pumping rate into 00°1 increases. 

Apparently this occurs because the populations of high-J rotational levels are 
greater at higher temperature. Consequently many pumping lines which are weak 
at low temperature produce more absorption as the temperature increases. The 
total flux emittted in the 9.4um and 10.4um bands is thus predicted to be 
temperature sensitive through the temperature sensitivity of the pumping rate. 

f. Future observations 


The fluxes produced by these models are compared to observed fluxes by 
Deming et al . (1983). The observed Mars flux is essentially 100% of the amount 
predicted by the models and the Venus flux is 74% of the amount predicted. The 
observed fluxes are a large fraction of the model prediction. We find this 
somewhat surprising since a large number of pumping lines were included in the 
calculations, and the various approximations (e.g. non-overlapping lin^s) should 
have the net effect of overestimating the emitted flux. Presuming that our 
calculations are largely correct, this means that the strength of this lOum 
emission is efficiently tied to the total instantaneous near-IR solar heating 
rate in the mesosphere. Nearly all of the bands important in near-IR heating 
(Dickinson 1972) are also important pumps for the lUnm emission. The strength 
of the lO^m emission should therefore be closely correlated with the 
instantaneous near-IR heating rate. Imaging of Mars and Venus, at high spatial 
resolution, and in a narrow frequency band containing the emission would be of 
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interest. If trace species such as water vapor contribute to this lOym emission 
we would expect such imaging to show non-uniformities. These variations would 
■ncctrssarily i)e related to variations Tn the radiative heating rate, since any 
process which contributes to the lOum emission will also leave CO 2 in Zv^, and 
this remaining energy is eventually thermal! zed. 

g) Speculations 

The laser amplification which occurs in these lines is not large by tne 
standards of laboratory laser technology. However, the volume of gas present in 
the region of level inversion is vast in comparison to laboratory samples, and 
the total solar pumping energy deposited in this region is certainly large in 
comparison to the power output of laboratory lasers. This suggests that it 
might be possible to utilize planetary mesospheres as a natural gain medium for 
the construction of large-volume solar-pumped lasers. In this regard, the fact 
that we calculate the long-path gain in these lines to be as large as 7% is 
significant. This gain is sufficient to overcome reflection losses, and thus 
.large orbiting mirrors could apparently be used’to define a lasing cavity in a - 
CO 2 pi anetary .mesosphere. We can use the present calculations to estimate the 
•power density that would be generated by such a cavity. In this case the rate 
equation (8) would be rewritten as a balance between the excitation of 00°1 
through the solar pumping term, and a dominant loss term by stimulated emission. 
In this limit the volume emission rate becomes equal to the solar pumping rate. 
In the case of Venus this corresponds to ‘^10 cm’ sec’ • The available path 
length is •/'2500 km at */‘13l3 km altitude. This gives an emergent power density of 
'^'2.5x10 lOym photons cm’ sec’ . This corresponds to only ‘^5x10’ watts cm’ , 
and so the power density obtained from such a cavity is not large. The reason 
for this is that the emission intensity is limited by the solar pumping, rate. 

Yet the solar pumping occurs through absorption in the Doppler cores of the 
near-IR CO 2 lines. Since these Doppler cores are narrow, the pumping rate is 
limited by the amount of frequency space available for absorption. Formation of 
a cavity with orbiting mirrors would therefore not increase the volume emission 
rate significantly over its present value, but it would serve to impose severe 
directionality and bandwidth-narrowing to this emission. The power, per unit 
solid angle, per unit bandwidth, emergent from such a cavity could achieve 
extremely large values. The possibility of producing very high specific 
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intensities from such large-scale lasers may have significance for the 
possibility of coirinumcation with extraterrestrial intelligence (CETI). Early 
technical studies ie.jg. Oliver 1974) roncJiided Jthat such coramunicacions ^were 
best conducted in the microwave region of the spectrum. Infrared communications 
at 10 Jim were considered, but found to be less favorable. If very high specific 
intensities can be obtained at 10 urn from large-scale mesospheric lasers, then 
these early studies should be re-assessed. 
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Figure Captions 


1 


Fig. 2 


Fig. 3 


Fig. 4 


Fig. 5 


^odel atmo3phBres adopted io ttie modelling of the IDum emission from 
the mesosphere of Mars. The models assume an isothermal mesosphere at 
120K and 170K. In situ temperature measurements by the Viking probes 
are shown in comparison. 

Model atmosphere adopted in the modelling of the lOMm emission from 
the mesosphere of Venus. The escape-to-space probability <c > for 

IP C 

the 4. Sum band of C 0^ is also shown. 

Schematic of important radiative and colli si onal processes in CO 2 , 
which are considered in our modelling of the lOum laser emission. The 
fL correspond to the notation used in the text for the viorational 
State populations (cm' ). 

Solar pumping rates into the (00°1) state of CO 2 for the 120K Mars 
model. The pumping by different frequency regions is shown 
separately, labelled by frequency in cm'^. The dashed line indicates 
the slope of the altitude dependence that would prevail for pumping by 
optically thin (unsaturated) absorption. 

Results of the model calculations for the 120K Mars model with the Sun 
at the zenith. The lower portion of the figure shows the calculated 
level populations (read left scale) for the upper (00°1) and lower 
(10°0) state for the lOum emission. The dashed lines give the 
population of these states under the assumption of local thermodynamic 

equilibrium. The shaded area shows tne region of level inversion. 

-3 -1 

Tne total number of photons emitted (cm sec , read rignt scale), in 
excess of thermal emission, in the 9.4 and 10.4um bands is shown as 
*^em' upper half of the figure shows the total optical depth in 
the core of the 10.33um R(8) line of ^^C^^02. This optical depth is 
calculated along a line perpendicular to the local zenith and is shown 
ras a function of the minimum altitude along this line. The dashed 
porrtion of the line indicates negative optical depths, i.e. 
amplification in the line. 
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Fig. 6 


Fig. 7 


Fig. 8 


Fig. 9 


Results of the model calculations for the 170K Mars model with the Sun 
at the zenith. See the caption of Fig. 6 for an explanation of the 
illustrated quantities. 

Results of the model calculations for the Venus model with the Sun at 
the zenith. See the caption of Fig. 6 for an explanation of the 
illustrated quantities. 

The modelled emergent intensity In the core of tne 10.33um R(8) line, 
shown as a function of angle for the 120K Mars model. The Sun is 
taken to be at the zenith, is defined as the angle between tne 
zenith and the direction of propagation, as me?<;iured at the top layer 
of the atmospheric model. 6p (top scale) is defined as the angle 
between the direction of propagation and the zenith, measured fran the 
peak of the emission at 67 km altitude. The intensity scale is 
normalized to unity at cos8g*1.0. The inset expands the relative 
intensity curve in the region of its peak. The dashed line shows the 
intensity obtained when laser amplification in the line is neglected. 

Total relative flux emitted in the 9.4 and 10.4um laser emission shown 
as a function of solar zenith distance, e^, for the 120K Mars model. 
The dashed line shows a cosine dependence. 
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TABLE I 

Principal near-infrared pumping 
given as '>^^24 v^r. 

bands. The vibrational 

BAND ORIGIN 
(an"') 

UPPER STATE 

LOWER STATE 

STRENGTH AT 
120K(cm molecule ) 

2336 

01111 

01101 

1.63x10“''^ 

2349 

00011 

OOOOl 

2.36x10’^^ 

3580 

11112 

01101 

1.70x10'^'’ 

3612 

10012 

00001 

2.55x10’^® 

3632 

10011 

00001 

3.95x20"®° 

3714 

10011 

00001 

3.70x10'^® 

3723 

mil 

01101 

2.41x10’®^ 

4978 

20012 

00001 

8.63x10"®° 

5100 

20011 

00001 

2.77x10"®° 

6972 

00031 

00001 

3.69x10“®^ 


state identifications are 

DJCrDENT SOLAR PLDX , 
AT 1.52 AU 

(Photons cm"^S“ ) 
1.19x10^^ 

1 . 62 x 10 ^^ 


2.01x10^^ 

2.23x10^^ 
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